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We study the effect of a transport current on the vortex-line lattice in isotropic type-II supercon- 
ductors in the presence of strong thermal fluctuations by means of 'driven-diffusion' Monte Carlo 
simulations of a discretized London theory with finite magnetic penetration depth. We calculate the 
current-voltage (I-V) characteristics for various temperatures, for transverse as well as longitudinal 
currents /. From these characteristics, we estimate the linear resistivities R xx = R yy and R zz and 
compare these with equilibrium results for the vortex-lattice structure factor and the helicity mod- 
uli. From this comparison a consistent picture arises, in which the melting of the flux-line lattice 
occurs in two stages for the system size considered. In the first stage of the melting, at a tempera- 
ture T m , the structure factor drops to zero and R xx becomes finite. For a higher temperature T z , 
the second stage takes place, in which the longitudinal superconducting coherence is lost, and R zz 
becomes finite as well. We compare our results with related recent numerical work and experiments 
on cuprate superconductors. 
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PACS numbers: 74.60.-w, 74.60. Ge 



I. INTRODUCTION 



The statistical mechanics of vortex lines in type-II su- 
perconductors has been the subject of intense study af- 
ter the discovery of the high-temperature superconduc- 
tors. Due to the large thermal fluctuations and pro- 
nounced anisotropies of these materials, thermal wan- 
dering of vortex lines leads to a rich phase diagram . 
Numerous experiments ]3|-[7j and computer simulations 
|js|-pTf were interpreted as evidence for the scenario of 
a vortex-lattice melting transition into a vortex-liquid 
phase. Most convincing evidence for the existence of 
a first-order phase transition separating the vortex lat- 
tice from a liquid phase comes from a recent calorimetric 
measurement of the specific heat in YBCO Q and from 
magnetization measurements in YBCO [|J as well as BiS- 
CCO ||. An explanation of the magnitude and temper- 
ature dependence of the observed characteristic entropy 
and magnetization jumps has been given in a recent work 
by Dodgson et al. Q . 

One interesting aspect of the melting transition is the 
question of whether it coincides with a complete loss of 
c-axis correlation (i.e. decoupling of the vortex lines into 
independent 'pancakes' in the Cu-0 layers). A recent 
transport measurement on untwinncd YBCO by Righi et 
al. p3[ indicates that just above the melting tempera- 
ture, the vortices are still correlated over a few microns, 
and become fully decoupled only at a distinctly higher 
temperature. 

On the theoretical side, the properties of the vortex 



liquid phase have been intensively studied. Some time 
ago Li and Teitel numerically studied the (uniformly 
frustrated) 3D XY model on a cubic lattice, and found 
a melting transition into a liquid with longitudinal su- 
perconducting coherence. The longitudinal supercon- 
ductivity, signalling a c-axis vortex correlation over the 
full system thickness, was found to be lost at a distinct 
temperature above the melting temperature. In later 
work including screening effects using the lattice Lon- 
don Model, a similar two-stage melting transition was 
found by Chen and Teitel Recently, a two-stage 

transition was also found by Ryu and Stroud Ji~8| using a 
3D XY model on a stacked triangular lattice, for a frus- 
tration lower than studied in earlier work by Hetzel et 
al. ||. When the intermediate liquid with longitudinal 
coherence would persist into the thermodynamic limit 
(very large system thickness), it would be the realization 
of the so-called line-liquid phase proposed in Ref. Q . 

In this paper we study the dynamical properties of 
the lattice London model and compare them with the 
behavior of equilibrium quantities as studied in Refs. 



ii. 



LATTICE LONDON MODEL AND MONTE 
CARLO METHOD 



The Hamiltonian of the isotropic lattice London model, 
at constant induction B, reads: 



1 



H = An 2 . J q^(Ri)q^(Rj)g(Ri - Rj) 



(i) 



(a) 



where J = <&Qii/(327r 3 A 2 ) (with A the magnetic pene- 
tration depth) and g(R) is the London interaction with 
Fourier components 



9{k) 



(d/xy 



(2) 



and k 2 = 



, KZ with k, ; = 2 — 2 cos fc„ 



rip = 0, 1, ...,Lfj, — 1). Here we assumed an L x X L y X 
L z lattice with periodic boundary conditions. At ev- 
ery dual lattice site R of our square lattice, the in- 
teger variable qu(R) denotes the vorticity or number 
of flux-line unit elements in the direction /_t = x,y,z. 
The q^{Ri) are subject to the continuity constraint 
X) e [<lti{Ri) - q^{Ri - e M )] = 0. Here Ri - e M runs over 
nearest neighbor sites of Ri. A is the magnetic penetra- 
tion depth and d the lattice constant. The Hamiltonian 
(Q) can be derived from the discrete version of the London 
free energy flTo|,[l4l 
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Monte Carlo sampling of the phase space for the vari- 
ables qfj,(R) at constant B is performed as follows. In 
the simulations, the initial configuration is prepared to 
contain the number of vortex lines corresponding to the 
value chosen for B. A Monte Carlo update step consists 
of adding at a given site a closed d x d square loop of unit 
vorticity with an orientation chosen randomly from the 
six possible ones. This scheme preserves the magnetic 
induction B with components B^ — -^>Ylj{q»(Rj))- 
(V = L x L y L z ). The standard Metropolis algorithm is 
employed to accept or reject the new configuration. For 
equilibrium calculations, one uses (|l|) to calculate the en- 
ergy change AE between the old and the new configura- 
tion. To simulate the presence of a uniform transport cur- 
rent j, we introduce an additional bias in the acceptance 
rates by subtracting or adding (AE)j = $ I 3 I d 2 jc 
to the energy change AE when an elementary loop is 
added whose normal vector is pointing in the same di- 
rection as j or —j respectively. We measure the magni- 
tude of the current in terms of the dimcnsionless quan- 
tity a = j&od 2 j (Jc). The dimensionless voltage V M is 
obtained by measuring the rate at which vortex jumps 
in the fj, direction are occurring. This 'driven-diffusion' 
method was used before in Refs. EG]. 



III. RESULTS 

Our calculations are carried out onal5xl5xl5 cubic 
lattice, with a magnetic field running in the z direction. 
We choose the filling fraction / = 1/15 and A = 5d, 
as in Refs. 
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_i|JT|]. Runs were taken consisting of 16,384 
(high currents) up to 262,144 (low currents) Monte Carlo 
sweeps through the lattice, half of which were used for 
equilibration. 
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FIG. 1. I-V characteristics, plotted as resistance versus 
current, (a) current along x direction, (b) current along z di- 
rection. Insets: I-V characteristics, plotted as voltage versus 
current, including the high current range. 

In Fig. 1(a) and (b) we show our results for the I-V 
characteristics at different temperatures, plotting the re- 
sistances V^/a M as a function of current for \x = x and 
H = z, respectively. In Fig. 1(a) the in-plane (/i = x) 
resistance is shown. For the lowest T shown (T = 1.5 
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in units of J/ks), the resistivity drops to zero below a 
finite current value (critical current) , and thus the linear 
resistivity 



R x 



V x 

lim — 

a.i^0 a x 



(3) 



is zero. For T — 2.0 no clear critical current is observed, 
but again we estimate R xx to be zero (or very small). 
For the higher T's studied, we find that the I-V char- 
acteristic is linear (constant resistance) for sufficiently 
small currents. For high currents, the I-V characteristics 
become independent of temperature (see insets in Figs. 
1(a) and (b)). In the high-current limit the voltage sat- 
urates at the value 1/6, corresponding to acceptance of 
all the Monte Carlo trial moves in the direction of the 
Lorentz force. 
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FIG. 2. Linear resistivities vs. T. Circles: in-plane. 
Squares: out-of-plane. Inset: temperature dependence of the 
vortex-lattice structure factor S = S(k s ), with k s the small- 
est reciprocal lattice vector of the ground state vortex lattice, 
and the helicity modulus y z , calculated at j = 0. 

In Fig. 2 we plot the linear resistivities R xx and R zz 
(estimated from the I-V characteristics of which repre- 
sentative ones are shown in Fig. 1) as a function of T. 
We distinguish three temperature regions with different 
dissipative properties. For low T, both R xx and R zz 
are zero. This is the flux-line lattice state, stabilized 
against small currents by the artificial pinning effect of 
the discrete mesh in our simulation. For intermediate 
T, R xx becomes finite but R zz is still zero. Finally, for 
high T, both transport coefficients become nonzero. For 
comparison, we include equilibrium (i.e. j = 0) results 
for the structure factor S and the helicity modulus 7 Z 
in the inset of Fig. 2. These equilibrium results were 
obtained before in Refs. The temperature at 

which the in-plane linear resistance becomes nonzero co- 



incides with the melting temperature T m of the flux-line 
lattice as found from the decay of the equilibrium struc- 
ture factor [ plf . Furthermore, the temperature at which 
the out-of-plane linear resistance R zz becomes nonzero 
coincides with the temperature T z at which the j = 
helicity modulus 7 2 vanishes (see inset). Here j z was 
calculated as explained in detail in Ref. |2q] . This is 
consistent with the interpretation of j z as a measure of 
longitudinal superconducting coherence or longitudinal 
superconductivity. Thus our transport calculations con- 
firm the two-stage melting picture found in equilibrium 
simulations of this model jl2 14|, in which there is an 
intermediate vortex-liquid regime with superconductiv- 
ity along the vortex lines. In this intermediate regime 
T m < T < T z the vortex lines are c-axis correlated over 
the full system thickness L z = 15 simulated. For T > T z 
the longitudinal correlation length becomes smaller than 
L z . At some even higher Td (decoupling temperature), 
the c-axis correlation should be lost completely. 



IV. DISCUSSION 



As our results are obtained for a finite system size, it is 
not clear from them whether the two-stage melting sce- 
nario persists into the thermodynamic limit. In fact, we 
believe that this will not be the case. In Refs. (rjlj it 
was found that the T z decreased with increasing system 
thickness L z . In Ref. jl9| , the L z dependence of T z was 
studied for the uniformly frustrated Villain Model. There 
it was also found that the intermediate temperature re- 
gion decreased in size with increasing L z . We interpret 
this finite-size behavior as a manifestation of a c-axis cor- 
relation length that, in a thick system, is finite and large 
just above the melting temperature and shrinks with in- 
creasing temperature until a decoupling temperature Td 
is reached. Thus, what these simulation results do pre- 
dict for real (i.e. thick) systems is that the vortex lattice 
melts into a liquid in which the vortex lines are c-axis 
correlated over many layer distances, i.e. far from de- 
coupled. Only for higher temperatures a crossover takes 
place into a decoupled regime above Td- This scenario 
agrees well with the transport experiments on YBCO by 
Righi et al. |p3[ , that indicate that just above the melt- 
ing temperature, the vortices are still correlated over a 
few microns ^?J, i.e. thousands of Cu-0 layer distances, 
and become fully decoupled only at a distinctly higher 
temperature. In contrast, measurements on BiSCCO by 
Doyle et al. were interpreted as evidence for a simul- 
taneous melting and decoupling. This is consistent with 
recent Monte-Carlo simulations of a strongly anisotropic 
modified 3D XY model by Koshelev [^0| , in which it was 
found that the longitudinal helicity modulus drops to 
zero close to the melting temperature already for L z = 40 
and larger. We note that such a behavior was also found 
for L z — 40 in a recent Monte Carlo study of an only 
slightly anisotropic 3D XY model, in which the melting 
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transition was investigated upon cooling 
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